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A quintessence scalar field in self-interacting Brans-Dicke theory is shown to give rise to a non- 
decelerated expansion of the present universe for open, flat and closed models. Along with providing 
a non-decelerating solution, it can potentially solve the flatness problem too. 
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^ I. INTRODUCTION 

a i 

The Standard Cosmological Model (SCM) can only describe decelerated universe models 
and so cannot reproduce the results coming from the recent type la supernovae observations 
upto about z ~ 1 [1] which favour an accelerated current universe. But as the SCM can give 
a satisfactory explanation to other observational properties of the present universe (e.g., 
primordial nucleosynthesis, extragalactic sources redshift, cosmic microwave radiation). The 
^ ■ recent extensive search for a matter field which can give rise to an accelerated expansion for 

the universe. This type of matter field is called 'quintessence matter' (or, shortly Q-matter) 
This Q-matter can behave like a cosmological constant [2] by combining positive energy 
density and negative pressure. So there must be this Q-matter either neglected or unknown 
responsible for this accelerated universe. At the present epoch, a lot of works has been 
done to solve this quintessence problem and most popular candidates for Q-matter has so 
far been a scalar field having a potential which generates a sufficient negative pressure. 
The quintessence proposal faces two types of problems [3] . One of these problems (referred 
as fine tuning problem, eliminate classically by Ratra and Peebles [2]), is the smallness of 
I ' the energy density compared to other typical particle physics scales. The other problem 

known as the cosmic coincidence problem is that although the missing energy density and 
matter density decrease at different rates as the universe expands, it appears that the 
initial condition has to be set so precisely that the two densities become comparable today. 
Quintessence has been proposed as that missing energy density component that along with 
the matter and baryonic density makes the density parameter equal to 1. A special form 
of quintessence field called the 'tracker field' has been proposed by Ratra and Peebles to 
tackle this problem [2] (see also ref. [4]). Since in a variety of inflationary models scalar 
fields have been used in describing the transition from the quasi-exponential expansion of 
the early universe to a power law expansion, it is natural to try to understand the present 
acceleration of the universe which has an exponential behaviour too, by constructing 
models where the matter responsible for such behaviour is also represented by a scalar 
field [5]. fnverse power law is the other potential [6] that has been studied extensively 
for quintessence models, particularly, for solving the cosmic coincidence problem. Re- 
cently, Bertolami and Martins [7] obtained an accelerated expansion for the universe in 
a modified Brans-Dicke (BD) theory by introducing a potential which is a function of 
BD scalar field itself. Very recently, Banerjee et al [8,9] also have shown that it is possi- 
ble to have an accelerated universe with BD-theory in Friedmann model without any matter. 

This paper investigates the possibility of obtaining a non-decelerating (q < 0) expansion 
for the universe in BD theory with scalar field which is minimally coupled to gravity and 
serves as the quintessence matter. 
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II. FIELD EQUATIONS AND SOLUTIONS 

The Brans-Dicke theory is given by the action [10] 



/ 



s= / V^g 



(j>R-^-<P ta <i>' a +L r , 



d A x (1) 



where <j) is the BD scalar field, uo is the dimensionless constant BD parameter and L m is 
the Lagrangian for all other matter fields. We have chosen the units 8ttG = c = 1. 

The matter content of the universe is composed of perfect fluid and a scalar field ip as 
the quintessence matter. We assume that the universe is homogeneous and we consider an 
anisotropic space-time with line-element 

ds 2 = -dt 2 + a 2 dx 2 + b 2 dfl 2 k (2) 
where a, b are functions of time only and 

dy 2 +dz 2 , when k — (Bianchi I model) 

dQ, 2 , = ^ d6 2 + sin 2 0d(/) 2 , when k — +1 (Kantowaski-Sachs model) 
d9 2 + sinh 2 9d(j) 2 , when k = — 1 (Bianchi III model) 

Here k is the spatial curvature index, so that the above three types [11] of models are 
Euclidean, closed and semi-closed respectively. 

Now the BD-field equations are 

b (bV k (p m + P J l (A 2 b4> 4> ( ) 



a b ab _ (p m +p i> ) 1 (A _M , J (4] 

a b ab 6 2 UJ \(b a b ) 6 6 U 



and the wave equation for the BD scalar field <p is 



^(^4j^=3T^[^-« + ^- 3 ^ (6) 

p m and p m are the density and the pressure of normal matter, and p^ are those due 
to the quintessence field given by 



P,=^ 2 + VW, P^ = \i> 2 -V{i>) (7) 



where V(ip) is the relevant potential. 

The wave equation for the quintessence scalar field ip is 
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•• a b\ ■ dV(tp) 

^(a +2 ^ = "^f (8) 
From the above field equations, we have the matter 'conservation' equation 



pm + \^ + 2 lj (Prn+Pm) =0 (9) 

Assuming that at the present epoch, the universe is filled with cold matter (dust) with 
negligible pressure, so using p m = 0, the conservation equation (9) gives 

Pm = £ (10) 

where p\ is an integration constant. 

Now we assume, the power law form of scale factors a, b and the BD scalar field <f> are 

a = ai t a , b^ht 13 , < j) = <j) 1 t 5 (11) 

where ai, bi, 4>i are positive constants and a, [3,5 are real constants with a + 2/3 > 3 
(for accelerating universe). 

From the field equations (3), (4) and (5) using (11) we have the expression for ip 2 as 

^ = _ _£l_ t - a -20 + fa [2/3 2 _ 2 a(a - 1) - u5 2 - 5(5 -I) -(a- 2^)5]^ (12) 

bf a\bf 

From (6) the potential V is given by 



v = ^L t S-2/3 _ * t -*-W + -M[{2u + 3)(a + 2(3 + 5-1)5 + 2{3 2 - 2a(a - 1) 

- w(5 2 - <J(5 - 1) - (a - 2[3)5}t s - 2 (13) 



The wave equation (8) for the quintessence scalar field ip can be written in the form 

After integration (14) the expression for V is 

fcfc (2a + 2/3 + 5) s _ w Pl w (2a + 4/3 + J -2) 

V ~ ~\ (5 - 2(3) * 2^&T 1 [ ' 

+ 5(5 - 1) + uj5 2 + (a- 2/3)5 - 2(3 2 } t 5 ~ 2 (15) 

Now, the consistency relations of the constants coming from the two identical equations 
(13) and (15) for the potential V are 



4a + 2(3 + 35 = 



(16) 
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(2w + 3)5(5 - 2) (a + 2/3 + (5 - 1) = (4a + 8/3 + 35 - 6)[2a(a - 1) + 5(5 - 1) 

+ w5 2 + (a - 2)5 - 2/3 2 ] (17) 

From these relations for an accelerating universe (a + 2(3 > 3), there are two possibilities: 
one in which k = and the second where k ^ 0. 

Case I : k = 

In this case , consistency conditions are (17) and 

a + 2/3 = c (18) 

where c is an constant (> 3). 

For solving (17) and (18), we may choose 5 = —2 and we have two possible solutions: 

(i) <* = /?=§, 5 = -2 

or 

(ii) a = 4 - c, /3 = -2 + c, 5 = -2 

For both solutions (i) and (ii) the expression for -0 2 (see eq.(12)) becomes 

iP 2 = — ^ t- c - 2{2oj + 3)0it~ 4 (19) 
This indicates that ui < —3/2 as ip 2 cannot be negative. 

For c = 4, one has 2|2w + 3|0i > and the deceleration parameter q = —j. 
In this case, equation (19) integrates to 

^ = ±j (20) 
where A 2 = — 2(2w + 3)<f>i — and the relation between V and ip becomes 

V = ViV 4 (21) 

where V\ being a constant, related to the constants a\,p\ etc. The model works for all 
time < t < oo with the condition 2|2w + 3\<fri > is satisfied. The value of uu is related 
to the other constants as follows: 



2(2u + 3)0! + ^ = 1±y/ w" V \ for solution (i) 



and 



2(2^ + 3)01 = ^-2^, (e = 0,l), for solution (ii) 

For c ^ 4, the model does not work for the whole range of time < t < oo. In this case, 
the deceleration parameter is q = - — 1. 



If c > 4, then the rate of acceleration is faster than q = —\ and from equation (19), 



ip 2 > restricts the validity of the model for t > t\ where 



2|2w+3|o7SJ^ 



(21a) 



Further if 3 < c < 4, then the universe will expands with an acceleration but with a rate 
less than q = —j and as before from equation (19) for real tp , the model works upto the 
time ti where 



*2 



Pi 



(216) 



For q = — i, the present age of the universe can be calculated from (5) as 



to 



„ „ A 2 V 1 A 4 ^ 1/2 
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[(H 2 ) + 2(H a ) (H b ) } 1/2 



where H a = ^ and Ht, = | . 
For large w limit, 



to 



\/^2ZJ 



1(^)0 + 2(^)0(^)0] 1/2 



where w is obviously a negative quantity. 
Case II : k ^ 0. 

In this case, consistency conditions are (16), (17) and 

a + 2/3 = 3 

Solving these three equations we have the following solutions: 
(i) a = 1, (3 = 1, 5 = — 2 for all values of w. 



(22) 



' 2uj+3 
2^-3 



provided for > 3/2 or w < —3/2. 

For the solution (i), the model works for a limited period of time, < t < t\ where 

" aibf 



h = 20i 



2 - (2w + 3) 



Pi 



(23) 



For an open universe i.e., for k = —1, (2ui + 3) < and \2u) + 3| > A. 
For a closed universe i.e., for k = +1, (2w + 3) < and \2u) + 3| < -k. 
For the solution (ii) the model works for the time, where t satisfies the equation 



2(2w + 3)0i 



it 



2u + 3 



2oj - 3 

provided for cj > 3/2 or ui < —3/2. 
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III. FLATNESS PROBLEMS AND ITS SOLUTIONS 

One important aspect of this model is that potentially it can solve the flatness problem. 
Now we make a conformal transformation [12] as 

g^v = <j> 9nv (25) 
In this section, we make the following transformations: 

dri=^/4>a, a=^/4>a, b = b, ip = In <f>, p m =<j)~ 2 p m , p v , = (f>~ 2 p , p , 

P m = <P~ 2 P m , = 4>^ 2 P^ (26) 
So the field equations (3) - (5) transformed to (after some manipulations) 

b) - 2 t- 2 jt + p = ^ +p * ) + —t 2 (j) (27) 



b'\ 2 a'V k _ (3 + 2lu) 



and 



I +2 uT + F=^ + ^ + —^^^ (28) 



W__aT_ k__a!V_ /6' x 2 
b a b 2 a b 



where '=ij- 

The BD scalar field in new version p is given by 
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We define the dimensionless density parameter £1 as 



Q = fi m +^+^ = J^ (31) 



(29) 



. = ^ffl 2 =, (30) 



where p — p m + p^ + p^ is the total density and fit are defined accordingly. 

Using (27)- (30), and combining the energy densities, we have the equation for the conser- 
vation for the total energy, 

p' + 3H{p + p) =0 (32) 

Here fl"=|^+2j^is the Hubble parameter in the Einstein frame and 7 is the net 
barotropic index defined as 

7^ = 7 m & m +% fi,+7* (33) 
From equations (28) and (32), we have the evolution equation for the density parameter 

as 
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Q' = 1)17^ +2(7- l)H b ] (34) 

where H = t and H .= \ . 

The individual 7^ 's are defined by the relation pi = (7, — l)pi- So the ratios ^- remain 
same in both frames. For our choices of matter, p m and p^ = p^ , so we have 7 m = 1 and 
7^ = 2. The other index 7^ is related by the equation 

% = P -^-^ = t^— (35) 

It has been shown that 7^ is varies with time. The equation (34) indicates that = 1 is 
a possible solution and this solution determines that (^^j < 0. This solution is stable for 
expanding universe (H > 0) if 7 < 3 with the relevant condition 

n m +4ft,<(2-37jft, (36) 

From the field equation (28), the curvature parameter ft k — —k/ b 2 vanishes for the 
solution = 1. So for BD-scalar field it is possible to have a stable solution corresponding 
to Cl = 1 and hence the flatness problem can be solved. 



IV. CONCLUDING REMARKS 



In this work, we have investigated the nature of the potential relevant to the power 
law expansion of the universe in a self-interacting Brans-Dickc (BD) cosmology with 
a perfect fluid distribution for anisotropic cosmological models. We have considered a 
non-gravitational quintessence scalar field ip with a potential V — V(ip). This scalar field 
in BD-theory is shown to give rise to an accelerated expansion for the present dust universe 
(where we have taken p m = 0) It is to be noted that at early stages of the evolution of the 
universe p m is non-zero. But if we take barotropic equation of state p m — (7 — l)p m , then 
equation (10) is modified to p m = ^ a p b l^ ■ 

We have presented a class of solutions describing non-decelerating universe for both flat 
(k = 0) and curved space-time (k 7^ 0). For k = 0, there are two possible solutions for 
different choice of the parameters. In both the solutions, the parameter lu must be negative 
(in fact 2oj + 3 is negative) to make the quintessence field real. The validity (on the time 
scale) of the solutions depends on the parameter c (defined in eq.(18)). For c = 4, the 
model works for all time < t < 00, while for c > 4, the model works for (t\,oc) where 
t\ is given by (21a). From this model we cannot predict the geometry of the universe 
before t\. Similarly for 3 < c < 4 the model is valid for t < t 2 i.e., from the begining to 
the time ti. Here also we cannot predict the nature of space-time of the universe after t = t 2 - 

In curved space-time (k ^ 0) there are also two possible solutions for different choice of the 
parameters. The time interval over which the solutions arc valid are given in equations (23) 
and (24). For close model (k > 0) the coupling parameter u is restricted by the inequality 
< 2uj + 3 < l/b\ while for open model w satisfies — 1/b 2 < 2uo + 3 < 0. From the above 
solutions, we note that the coupling parameter lu may be positive or negative (with some 
restrictions). Hence for all solutions the choice of to is not in agreement with the observations. 

Finally, for non-decelerating solution it can also potentially solve the flatness problem 
(without any restriction on the parameters) and it has been shown that Q = 1 could be a 
stable solution in this model. 
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